function [dataXrtm,erroX]=Sqretiming1(dataXrt)
% Resample 4Sa/s

dataXrt=reshape(dataXrt,1,[]);
N=4;
L=length(dataXrt)/2;
dataXrt=resample(dataXrt,N,2);
%sqrt
dataXk=abs(dataXrt).^2;
% figure(301);pwelch(dataXk);
indexp=(0:L*N-1)/N;

indexp=exp(-1i*2*pi*indexp);

dataXp=dataXk.*indexp;
% erro
Xp=sum(dataXp);
erroX=-1/(2*pi)*angle(Xp);
Oindex=1:N*L;
Rindex=1:2:N*L;

XRindex=Rindex+erroX*N;
%retiming

IX=interp1(Oindex,real(dataXrt),XRindex,'cubic');
QX=interp1(Oindex,imag(dataXrt),XRindex,'cubic');
IX(isnan(IX)) = 0;
QX(isnan(QX)) = 0;
dataXrtm=complex(IX,QX);




